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Axisymmetric solitonic states {ehiral skyrmions) have been predicted theoretically more than two 
decades ago. However, until recently they have been observed in a form of skyrmionic condensates 
(hexagonal lattices and other mesophases). In this paper we report experimental and theoretical 
investigations of isolated chiral skyrmions discovered in PdFe/Ir(lll) bilayers two years ago (Science 
341 , 636 (2013)). The results of spin-polarized scanning tunneling microscopy analyzed within 
the continuum and discrete models provide a consistent description of isolated skyrmions in thin 
layers. The existence region of chiral skyrmions is restricted by strip-out instabilities at low fields 
and a collapse at high fields. We demonstrate that the same equations describe axisymmetric 
localized states in all condensed matter systems with broken mirror symmetry, and thus our hndings 
establish basic properties of isolated skyrmions common for chiral liquid crystals, different classes 
of noncentrosymmetric magnets, ferroelectrics, and multiferroics. 

PACS numbers: 12.39.Dc; 68.37.Ef; 75.70.Ak; 75.70.-i 


I. INTRODUCTION 


Long-period homochiral magnetization modulations 
{helical phases)^ and axisymmetric solitonic patterns 
{vortices or skyrmions'P^ are two types of unconven¬ 
tional magnetic states attributed solely to magnetic 
compounds with broken inversion symmetry and distin¬ 
guish them from common (achiral) magnetic materials 
(Figs. Both, extended chiral modulated phases 

and localized skyrmionic states are stabilized by specific 
Dzyaloshinskii-Moriya (DM) interactions arising in chiral 
magnets owing to their crystallographic handednes^. In 
the micromagnetic energy functionals of noncentrosym¬ 
metric ferromagnets these interactions are described by 
energy contributions linear in the first spatial derivatives 
of the magnetization M {Lifshitz invariants)^ 
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Axisymmetric localized structures (Fig. are re¬ 
lated to multidimensional topological solitons with non¬ 
singular internal structure and finite energ}!^. These 
particle-like objects are of special interest in fundamen¬ 
tal physics and mathematics^^. In most nonlinear phys¬ 
ical systems, multidimensional solitons (skyrmions) can 
exist only as dynamic excitations while static solutions 
are unstable and collapse spontaneously into topological 
singularities^. 

In nonlinear field theory, the existence and stability of 
skyrmion solutions is provided by special terms in the 
energy functionals. More than five decades ago T. H. 
Skyrme introduced into the nonlinear field model an in¬ 
teraction term with higher order spatial derivatives that 
stabilize two- and three-dimensional topological nonsin- 




FIG. 1. (color online). Axisymmetric isolated skyrmions: 
(a) in cubic helimagnets and uniaxial ferromagnets with Dn 
symmetry; (b) in uniaxial ferromagnets with Cnv symmetr/^. 


gular solitons (now commonly addressed as skyrmionsj^. 
Since that time, field theorists have been intensively in¬ 
vestigating this family of solitons (sk yrmion s) within the 
Faddeev-Skyrme and kindred models^^^^^. 

Lifshitz invariants of type Q provide the only known 
alternative to the Skyrme mechanism that y ield regular 
solutions for axisymmetric skyrmions^E^ESl. These in¬ 
variants arising in noncentrosymmetric condensed matter 
systems (including chiral magnets, liquid crystals, multi¬ 
ferroics, and nanolayers of magnetic metals with interface 
induced Dzyaloshinskii-Moriya interactions) introduce a 
unique class of materials where mesoscopic skyrmions can 
be induced and manipulated. 

In a broad range of applied magnetic fields and tem- 
peratu res iso lated skyrmions condense into hexagonal 
latticej^^^^ ^ or other types of two-dimensional modu¬ 
lated state^^^^. During the last years, intensive ex¬ 
perimental efforts have been undertaken to find indica¬ 
tions of hexagonal skyrmion lattices in different groups 
of chiral ferromagnets (see and bibliography 

iiP). Particularly, direct observations of skyrmion lat- 
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tices have been reported in free standing nanolayers of 
cubic helimagnets irP^ (and the following papers of this 
gro^ppmU^ ^ These results reveal axial symmetry and ho- 
mochirality of the embedded skyrmions, and observed 
properties of skyrmion lattices were found to be in close 
correspondence with theoretical results. To date the 
LTEM studies of confined cubic helimagnets have focused 
on the skyrmion condensates (skyrmion lattices and 
clusters p^ l ^^ l ^^ l ^^ ^. Spin-polarized scanning tunneling mi¬ 
croscopy (SP-STM) has been able to identify isolated 
skyrmions in the saturated states of PdFe/Ir(lll) filmP^, 
and subsequently resolve their internal structur^^. 

In this paper we present detailed experimental 
and theoretical investigations of axisymmetric isolated 
skyrmions in thin magnetic films. 

In the theoretical part we develop a consistent theory 
of chiral skyrmions in thin magnetic layers (Sec. IIA). 
In Sec. IIB we apply the qualitative theory of differen¬ 
tial equations to expound main features of isolated chiral 
skyrmions and elucidate their physical nature, investigate 
the conditions of the elliptical instability at low fields and 
calculate within the discrete model the skyrmion collapse 
field. In Sec. IIC we construct the phase diagram of the 
solutions for isolated skyrmions. 

In the experimental part we present the detailed evolu¬ 
tion of isolated skyrmions in PdFe/Ir(lll) bilayers from 
the strip-out at low fields to the collapse at high fields. 


II. THEORY 



FIG. 2. (color online). Basic modulated phases in chiral fer- 
romagnets: one-dimensional helicoids (a) and cycloids (b) and 
two-dimensional skyrmion lattices (c,d). Bloch-type modula¬ 
tions (a,c) arise in cubic helimagnets and ferromagnets with 
Dn symmetries; Neel-type modulations (b,d) are attributed 
to uniaxial ferromagnets with Cnv symmetriefl 


fornffl 


w = A(grad m)^ + WD{m) — K{m ■ n)^ 

— ' n —-jaoMm ' (2) 

where A is the exchange stiffness constant, K is the 
uniaxial anisotropy constant, is the demagnetizing 
field. 


A phenomenological theory of chiral modulations in 
noncentrosymmetric magnetic crystals has been devel¬ 
oped by I. Dzyaloshinskii in 1964^. These papers also 
include analytical solutions for one-dimensional chiral 
modulations {helicoids and cycloids). Theory of isolated 
skymions and skyrmion lattices in bulk no ncen trosym- 
metric ferromagnets has been developed irP^. Theo¬ 
retical investigations of chiral modulations in bulk and 
confined noncentrosymmetric ferro- and antiferromag- 
nets have been carried out in many of the papers dis¬ 
cussed in RefP. 


m = M/|M| = (sin ^ COS 7/^, sin ^ sin t/;, COS (3) 

is the reduced magnetization, n is the unity vector di¬ 
rected perpendicular to the film surface. 

The Dzyaloshinskii-Moriya energy density wd is com¬ 
posed of Lifshitz invariants 0 : 
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A. The micromagnetics of chiral modulations 

1. Energy functional and symmetry 

In this paper we investigate isolated skyrmions in a 
thin layer of a noncentrosymmetric ferromagnet. As a 
model we consider a thin plate infinite along the x— and 
y— axes and of thickness L along the z— axis. In the 
following sections we specify the model and discuss its 
limitations. For a film of a noncentrosymmetric uniaxial 
ferromagnet in the applied magnetic field perpendic¬ 
ular to the film surface, the micromagnetic energy den¬ 
sity written within terms quadratic in the components of 
the magnetization vector M has the following standard 


The functional forms of energy density wd are deter¬ 
mined by crystallographic symmetry of a noncentrosym- 
metr ic magnetic crystal and are listed in Eqs. ( |A.1[ ), 
(A.2). Lifshitz invariants 0 javour spatial modulations 
with a fixed rotation sense along the Xk direction^. A 
competition between the chiral energy wd and other en¬ 
ergy contributions leads to the formation of isolated chi¬ 
ral stateJ^ and spatially modulated magnetic phaseJ^^^. 

The Euler equations for energy functional 0 together 
with Maxwell’s equations. 


rotH('^) = 0, div = 0, 


( 5 ) 


yield solutions for different types of chiral modulations 

(Figs, m Bp). 
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2. Demagnetization effeets 

Generally the equilibrium modulated patterns m(r) in 
a chiral magnet are derived by numerically solving the 
above set of nonlinear different ial e quations including 
non-local stray field calculation^^^. Contrary to soft 
magnetic materials where demagnetization fields suffi¬ 
ciently influence the equilibrium magnetic stateJ^, in 
chiral magnetic materials the DM interactions strongly 
suppress these effect^. As a result in many practi¬ 
cal cases a magnet ostatic problem is reduced to analyti¬ 
cal solutionP^^^^^, and the stray-field energy can be ex¬ 
pressed as local energy contributions in energy functional 

It was also found that for one-dimensional modulations 
and two-dimensional axisymmetric struct ures, the inter¬ 
nal stray-field energy has a local character^^. Particu¬ 
larly, for ferromagnets with Cnv symmetry the internal 
stray-field energy can be taken into account by the fol¬ 
lowing redefinition of the anisotropy constant. 


K^K + Kd, Kd = ^ioMy2. 


3. The equations for axisymmetie skyrmions 


(6) 


We introduce cylindrical coordinates for the spatial 
variable r = (r cos (p, r sin cp, z) and consider magnetic 
patterns homogeneous along the 2 ;-axis with the magneti¬ 
zation antiparallel to the applied field in the center {0 = n 
for r = 0) and approaching the parallel orientation when 
the distance from the center approaches infinity ^ 0 
for r ^ oc). For 0{p^ip), the energy functional 

^ is reduced to the following form: 


w = A 


^ f V’r + 


■Wd 


- K cos^ 0 - COS 0 — poMm • H 


(d) 


( 7 ) 


and the Dzyaloshinskii-Moriya energy functionals 
WD{0,2p,r,ip) are listed in Eqs. ( A.1^(|A.2|). 


The equations minimizing energy 
ally symmetric solutions. 


A^( 

gy 0 


include rotation- 


0 = 0(p}, = H(‘') =H(‘')(p). 


(8) 


Analytical solutions for uniaxial noncen- 

trosymmetric ferromagnetand cubic helimagnets (Figs. 
0[1) are listed in Eq. dVTt - 

To date, only two types of skyrmionic states from this 
list have been identified in chiral ferromagnets by di¬ 
rect experimental observations: skyrmionic patterns with 
Bloch-type modulations (Eig. [^a) 

m = sin 0{p) -h cos 0{p) (9) 

have been observed in free standing nanolayers of cubic 
helimagnets (see and skyrmion lattices with 

Neel-type modulations (Eig. 0 b) 


have bee n observe d in Ee/Ir(lll) and PdEe/Ir(lll) 
nanolayer and in the rhombohedral ferromag- 

net GaV 408 with Gs^ symmetr^EU. 

The first direct observations of isolated skyrmions have 
been reported in PdEe/Ir(lll) nanolayers^. These chi¬ 
ral solitonic structure s hav e been investigated in a broad 
range of applied field^^^^. 

After integration with respect to cp, the total energy 
for an isolated skyrmion of Bloch- and Neel-type in an 
applied magnetic field perpendicular to the film surface 
can be reduced to the following form: 


r =2Tr fP f{0, r)rdr. 


( 11 ) 


Here f{0,r) = w{0,r) — re(0) is the difference between 
the skyrmion energy density and that of the saturated 
state, re(0) = —K — poMH : 

f{0, r) = A ^ sin^ 0^ — D + - sin 6> cos 0 

+ K sin^ 0 -1- pqH (1 — cos 0) . (12) 

In Eq. H = Hz is the perpendicular component of 
the internal magnetic field that differs from the applied 
external field due demagnetization field of the film 

surfac^^. Eor rather thick films {d > Vg) H = — 

PqM and for ultrathin films H ^ 

The Euler equation for energy functional (p!^. 


(13) 


A (Orr + \ siu 0 COS | — siu^ 0 

\ r Jr 

— K sin 0 cos 0 — poMH sin 0 = 0, 
with boundary conditions 

6>(0) = TT, 0{oo) = 0, 


(14) 


yields the equilibrium structure of isolated axisymmet¬ 
ric skyrmion^^. Note that for Neel-type skyrmions K 
includes the stray energy contribution 
Dimensionless variables 

p = 27rrlLD, h = H/HD, k = K/K^, (15) 

are commonly used in recent papers to describe mod¬ 
ulated states in uni axial chira l ferromagnets and cubic 
helimagnets (see e.g.ESlMEniMl)^ Here we use the c harac - 
teristic parameters of a uniaxial chiral ferromagnefPESl 

47rA 

= = 7^0 = ^. (16) 

Lb is the period of a helix at zero field and zero 
anisotropy, Hb is the saturated field and Kq is the criti¬ 
cal anisotropy ( A.13|). 

With variables (15), the equation for axisymmetric 
skyrmions (13) is reduced to the following form: 


^ 1 . ^ ^ 2sm^6> 

Opp H— ^ sm 6 cos 0 H- 

p p^ p 

— k sin 0 cos 0 — h sin ^ = 0, 


m = e p sin 0{p) -h cos 0{p) 


(10) with boundary conditions (14). 


(17) 
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FIG. 3. (color online), (a) Typical localized solutions of the boundary value problem with magnetization profile 0{p). The 
first three excitation modes Ci{p) with positive eigenvalues (b) for the solution for KjK q = 2.8. “Shooting trajectories 0{p) 
of the Cauchy problem (c) and corresponding phase trajectories 0p{0) (d). 


B. Solutions for axisymmetric skyrmions 


The equilibrium skyrmion profiles 0{p ) are derived by 
solving the boundary value problem ( p^ and (14) with a 
finite-difference methocP. Typical solutions of Eq. (13) 
are plotted in Fig. and the existence areas for isolated 
skyrmions are indicated in the phase diagram of the so¬ 
lutions (Fig. [^. 

The solutions 0{p) are linear near the skyrmion axis 
{{tt — 0) oc p for p ^ 1) and decay exponentially 
at high distances from the center {p 1) 0 oc 
exp {-pVk + h) /yp. 

Usually the functions 0{p) have arrow-like shape with 
the steepest slope at the center of the skyrmion (r = 
0). They transform into bell-shape profiles only near the 
critical line Hei . In micromagnetism, the characteristic 
size of a localized magnetization profile 0{p) is defined 
oJMl 


rs=ro- Oo {d9/dr)^l^^ , (18) 

where (ro,^o) is the inflection point of the profile 0{r) 

(Fig. I^a). 

The basic properties o f the solutions for Eq. 
have been investigated irP^. Theories of static chi¬ 
ral skyrmions in different classes of bulk and confined 
chiral magn ets h ave been developed in a number of 
studies (e.g.j ^^ l ^^ KHj^ Numerical solutions for isolated 
skyrmions in nanodots and other confined chiral mag¬ 
nets have been derived in a large number of recent 
works (e.g.lilBSHill)^ Also, dynamical properties includ¬ 
ing current-induced movement of skyrmions have been 
intensively investigated by numerical simula tions of the 
Landau-Lifshitz-Gilbert equation (e.g., Refs.SS^^^^^ and 
the bibliography in a review paper^. The results of 
these numerical simulations demonstrate a rich spectrum 
of magnetic states characteristic for chiral skyrmions and 
various scenari os of th eir evolution under the influence 
of applied fleld d^^ * ^^ ^. Particularly, in confined uniaxial 


helimagnets the applied field induces modulated textures 
with different number of skyrmions, elongated, and half 
skyrmion^^. 

The results of numerical simulations for stationary 
and moving skyrmions, however, still require substan¬ 
tial analytical analysis and physical comprehension. The 
qualitative theory of nonlinear differential equations to¬ 
gether with other analytical methods provide effective 
tools to gain important insight into the physics of chiral 
skyrmions and establish mathematical relations between 
them and other types of magnetic solitons. 


1. Visualization of solutions on the (0,0r) phase plane 


Solutions 0{r) of the boundary value problem (14) can 
be derived by solving the auxiliary Cauchy initial value 
problem for equation (HI, 


^(0) = TT, 6>r(0) = —a. 


(19) 


For illustration we consider the Cauchy problem given by 
(13) and (19) for iT = 0 and >c = 'kDI{A\/AK) = 0.8 
( A.16[ ). The calculated profiles 0{r^a) and the corre¬ 
sponding curves 6>r(6>) in the interval [0.4 < a < 4.0] 
are plotted in Fig. (c),(d). Most of curves 0{r,a) 
oscillate near lines 0 i ^2 = =k7r/2, the maximum values 
of Wa = and the corresponding profiles 0r{0) 

spiral around the attractors^ points (±7r/2,0). Among 
these curves there is a singular line (with a = 1.62471) 
which ends in the saddle point (0, 0) and, thus, repre¬ 
sents a solution of the boundary value problem for iso¬ 
lated skyrmions. 

The visual representation of the solutions for the aux¬ 
iliary Cauchy problem ( p!^ , (19) as parametrized profiles 
6>(r, a) (Fig. [^(c)) and 0^0) curves in (6>, Or) phase plane 
(Fig. I^(ci)) reveal mathematical regularities in the for¬ 
mation of the localized states. 
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FIG. 4. (color online). In the phase diagram in variables k 
and h the existence area of metastable isolated skyrmions is 
restricted by the strip-out critical line hei(k). The inset shows 
the regions of global stability of the modulated (helicoidal and 
skyrmion lattice) and the spatially homogeneous saturated 
phases (for details sej^. 


To demonstrate a crucial role of the DM interactions 
in the stabilization of chiral skyrmions, in the following 
we compare the phase portrait in Fig. [^(d) with special 
cases of model ( 12 ) with D = 0 . 


Isotropic ferromagnets [D = K = H = 0). The Eu¬ 
ler equations for energy functional of an isotropic ferro- 
magnet w = 74 (grad m)^ yield rigorous analytical solu¬ 
tions for axisymmetric skyrmions ^(r), derived by 

Belavin and Polyako\0^ 


^ = Nip-\-a^ tan (^/ 2 ) = (J/r)^, ( 20 ) 



FIG. 5. (color online). Magnetization profiles for Belavin- 
Polyakov instantons (a) and the corresponding phase por¬ 
traits of the solutions (b). 


2. Derrick scaling identities and a virial theorem for chiral 
skyrmions 


Analysis of skyrmion energy T 0 under scaling 
transformations offers further important insight into the 
physics of chiral skyrmions. We consider a family of func¬ 
tions 'd(r) = 'dirjrf) obeying the boundary conditions 
( p^ . Here 77 > 0 is an arbitrary constant describing uni¬ 
form compressions (0 < 77 < 1 ) or expansions (77 > 1 ) of 
profile 'd(r). For rescaled functions 'd(r) = 'd(r/77), the 
skyrmion energy T ( 11 ) can be expressed as a function 
of 77 : 


J^{r])=£e-£D'n + £o'>f- (22) 

The values of the exchange {£e)i Dzyaloshinskii-Moriya 
{£d)^ and potential (Dq) energy contributions for profile 
'd(r) 0 are given as follows: 

fe = 27rA ^ sin^ = Aa\^ 


where a and > 0 are arbitrary values and N are positive 
integers. The energy 0 for solutions ( [^ tFo = Sir AN, 
does not depend on values S and o®. For = 1 a set 
of magnetization profiles 0{r/6) ( 20 ) and phase portrait 
trajectories 0r{0) 


0 = 2arctan((5/r), SOr = —2 sin^( 6 >/ 2 ), ( 21 ) 

are plotted in Fig. For (5 > 0, the curves 0r{0) start 
in points (tt, —2/(5) and end in the saddle point ( 0 , 0 ). 
However, any anisotropy or magnetic field will destabilize 
this solution. 

Uniaxial centra symmetric ferromagnets {D = 0). In 
this case Eq. 0 has no stable solutions for isolated 
skyrmions. Eor H > 0 all phase trajectories OJO) spiral 
around attractor (— 7 r/ 2 , 0 ). Eor H <0 Eq. ( p^ has ra¬ 
dially unstable solutions for isolated skyrmions as proved 
by Derrick-Hobart theorem (Eor details sed^^^. 


£d = 27r|D| + I sin 79 cos 7 ^^ = \D\a 2 , 


£0 = 27r [K sin^ 79 + fioMH (1 - cos 79 )] ^d^ (23) 

or fo = Kas ^ poMHa^, where ai are the numerical 
coefficients given by the values of the integrals in Eqs. 

([^. 

Eq. ([^ shows that the DM energy plays a crucial 


role in stabilizing skyrmion^^^. In centrosymmetric fer¬ 
romagnets {£d = 0 ) isolated skyrmions are unstable with 
respect to compression and collapse into a singular line 
(77 ^ 0) {Derrick-Hobard theoreuP). Skyrmion solutions 
that minimize the free energy ( [ 2 ^ only occur for nonzero 
Dzyaloshinskii-Moriya energy contributions. 

Ansatz solutions. Potential fF{r]) (22) has a convenient 
form for analysis of skyrmion solutions with trial func¬ 
tions of type 79 = 79 ( 77 / 77 ) that obey the boundary condi¬ 
tions (p!^. Particularly, a linear ansatz 


79 = 7 r[l — (r/ 77 )] (r < 77 ), 79 = 0 (r > 77 ), (24) 
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has been used in RefP to introduce the phenomenon of 
chiral skyrmions. The ansatz, 


d{r/r]) = 4 arctan [exp (—r/ 77 )] , 


(25) 


based on solutions for isolated 360° Bloch wallP^ pro¬ 
vides a good fit to the solutions of Eq. (13). In Ref.l^, 


magnetization profiles for isolated skyrmions have been 
fitted by a combination of functions of type (25). For the 


trial function ^{r/r]) in Eq. ( [2^ , the total energy (22) 
can be written as 

J^{r])/{27r) = 4.31A^{l.b9K^1.39iJ.oMH)r]‘^-3mDr]. 
For zero anisotropy {k = 0) this ansatz yields the transi¬ 
tion field into the skyrmion lattice hg = Hg/Ho = 0.760 
(cf. with the rigorous value hg = 0.801 and hg = 0.675 
for the linear ansatz (24lP). 

For J-{r]) (22) the equilibrium skyrmion size is 


L 


Vo = 


D 


0^2 


27r ask2a4h^ 


(26) 


expressed as a ratio of the Dzyaloshinskii-Moriya to the 
potential energy contributions for the trial function. 

The virial theorem for isolated axisymmetric 
skyrmions is derived by integration of the Euler 
equation (13). Partial integration leads to the following 


virial relation between the equilibrium values of the 
potential and DM energieJ^ Eq = 2£d where fo, Ed 
are the integrals (23) calculated for the solutions of Eq. 


3. Radial stability and collapse at high field (discrete model) 

The stability of the solutions 0{p) of the boundary 
value problem ( p^ , ( p^ under small radial distortions 
i{p) (^(0) = = ^has been investigated iiP. This 

problem is reduced to the ^ectral problem for the per¬ 
turbation energy functionaP. By numerically solving the 
eigenvalue problem for this functional, the radial stabil¬ 
ity of isolated skyrmions has been established in a broad 
range of the control parameters (/c, Contrary to 

magnetic bubbles, which collapse with finite radii at cer¬ 
tain critical fields^, the solutions of the boundary value 
problem © and ( p^ (Fig. exist at arbitrary high 
fields. In increasing fields their sizes gradually decrease 
and asymptotically approach zero. 

The continuum model ([^, however, becomes invalid 
for localized solutions with sizes of few lattice con¬ 
stants. In this region we investigate solutions for chi¬ 
ral skyrmions within the discrete models. We consider 
classical spins, S^, of unit length on a two-dimensional 
square lattice with the following energy functionai^ E = 
Eq E]j where 

Eo = -J ■ “)'](27) 

<i,j> i 

and the Dzyaloshinskii-Moriya energy equals 




FIG. 6. (color online). Micromagnetic energies of an iso¬ 
lated skyrmion (a) and a bubble domain (b) as functi on of 
their sizes for selected values of the applied magnetic fieldP^. 
Isolated bubbles collapse at critical field Hbc with finite radius 
rhc- Isolated chiral skyrmions exist at very high fields with¬ 
out collapse. The equilibrium skyrmion sizes p as funct ions 
of the applied field calculated for different values of po (30) 
indicate the collapse of chiral skyrmions (c) (p is defined here 
as a diameter of a circle encompassing a skyrmion core area 
with ruz < 0.995). 


Ed = -D y] (Si X Si+* • X + Si X Si+j • y) (28) 

for Bloch-type modulations, and 

Ed = -D y](Si X Si+s - y-SiX Si+y ■ x) (29) 


for Neel-type modulations (< > denotes pairs of 

nearest-neighbor spins). 

For a helix = (cos^^, sin^^, 0) propagating along the 
x-axis at field and anisotropy (H. = K = 0), model (27) 
is reduced to 


E 


= E - 7)sin(6»i - 9i+$)], (30) 


and yields the equilibrium period po = 27r/ arctan(D/J) 
(po is the number of magnetic ions corresponding to AO = 
27r). 

The calculated equilibrium skyrmion diameter p as a 
function of the applied field indicates the collapse of the 
skyrmion core at certain finite fields hc{po) (Fig. §(a)). 
The critical field hc{po) increases without limit with in¬ 
creasing Po (Fig. (b)) and, thus, signifies a transition 
from the discrete model to the continuous model. 


4- Elliptic instability (strip-out) at low fields 

Isolated skyrmions exist as metastable states above the 
critical field hg{k) (Fig. |^(a), (b)). Below this line the 
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FIG. 7. (color online). Collapse (He) and strip-out or elliptic 
instabil ity (Hei) critical fields calculated within th e di screte 
model ^7\ for /c = 0 and different values of po (30). At 
the dashed line Hs{po) the isolated skyrmion energy equals 
zero, and below this line skyrmions can condense into the 
hexagonal lattice (a). Calculated distributions of the magne¬ 
tization of the skyrmion core for k = 1.315, po = 24 and for 
h = 0.42 (b) and h = 1.27 (c). 


energy {11) becomes negative and skyrmions tend to 
condense into a hexagonal lattic^. However, if the for¬ 
mation of skyrmion lattices is suppressed (as in PdFe/Ir 
(111) films^ isolated skyrmions continue to exist be¬ 
low the critical line hs{k) (with the skyrmion core en¬ 
ergy density lower than that of the surrounding satu¬ 
rated state). At the same time isolated skyrmions have a 
tendency to elongate and expand into a band with heli- 
coidal or cycloidal modulations and eventually to fill the 
whole space, since the spiral state represents the mini¬ 
mum with lower energy as compared to the local minima 
with the metastable isolated skyrmions. These {ellip¬ 
tic) instabilities are similar to ” strip-out” instabilities of 
isolated magnetic bubbles at a certain critical fielcP^ ob¬ 
served in common “bubble-domain” film^^ and in mag¬ 
netic nanolayers with perpendicular anisotrop^P^. For 
chiral skyrmions, the elliptic instability fields are cal¬ 
culated from the stability analysis of the skyrmion energy 
(11) with respect to (elliptic) perturbations of typ^^ 


p = p-\-sr]{p)cos2(p, = i;-\-C{p,ip), (31) 


{e ^ 1). For isolated Bloch-type skyrmions the calcu¬ 
lated critical line hei{k) { 0 < k < ka ) is plotted in Fig. 

These results are close to earlier calculations for stray- 
field free elliptical distortions (31) with ansatz functions 


r]{p) = sinO /(1 + asin^) optimized with respect to the 
parameter 


Within the discrete model (27) the critical field hei 


has been calculated for zero anisotropy (/c = 0) and for 
6 < Po < 30. Fig. [^(a) shows that the strip-out field h^i 
essentially decreases with the decreased size of skyrmions 
what can be beneficial for possible application of such 
skyrmions. However, the existence region of these iso¬ 
lated skyrmions is restricted by the lower field of collapse 

hn. 



FIG. 8. (color online). Magnetization profile 0{r) for an 
isolated skyrmion (1) (Inset) derived from SP-STM data for 
the applied field — 1.11 T. The solid line is the solution 
of Eq. ( pTf for k — 1.315 and h = 0.321, Vg is the skyrmion 
core radius defined by Eq. (18). 


5. The k - h phase diagrams 


In this section we consider the existence area for iso¬ 
lated skyrmions in the magnetic phase diagram (Fig. [^. 
The energy functional for uniaxial chiral ferromagnets 
(© depends on the two independent control parameters, 
the reduced values of the applied field, h and uniaxial 
anisotropy, k (15). The magnetic phase diagram in vari¬ 
ables k and h collects all possible solutions for model (§. 
The calculated phase diagram in the inset of Fig. shows 
the existence areas of the cycloids and skyrmion lattices 
and the transition lines between these modulated phases 
and the saturated state. The phase diagram indicates 
the critical fields at zero anisotropy, the bicritical point 
B (1.90, 0.10), and the critical point A (2.67, OP (for 
a detailed description of this phase diagram see Ref.tl^. 
Fig. shows critical lines for isolated skyrmions (results 
of the continuum model <§)• Isolated skyrmions con¬ 
dense into a skyrmion lattice when the applied magnetic 
field decreases to the critical value hs{k). However, iso¬ 
lated skyrmions can exist as localized objects below the 
critical line hs{k) and strip-out into helicoids at the crit¬ 
ical line hei. 


III. EXPERIMENT: ISOLATED SKYRMIONS 
IN PDFE/IR(111) NANOLAYERS 

Sample preparation and spin-polarized (SP)- STM ex¬ 
periments were performed in a multi-chamber UHV sys¬ 
tem at a base pressure of 5e-ll mbar. Details of the 








































FIG. 9. (color online). The extension of elliptic deformations in isolated skyrmions in decreasing applied magnetic fields (a-f) 
is terminated by the formation of cycloid patterns (g,h). 


sample preparation can be found in Ref.^^. We use an¬ 
tiferromagnetic bulk Cr tips to minimize magnetostatic 
interactions between tip and sample. The SP-STM mea¬ 
surements were performed at T=4.2 K in perpendicular 
magnetic fields of -3 to +3 Tesla. We repeatedly scanned 
the same sample area while continuously sweeping the 
magnetic field at a speed of 12.8 mT/min, resulting in 
a series of images with a field difference of AB=S7 mT. 
Constant current images and maps of differential conduc¬ 
tance (dl/dU) were measured simultaneously by a lock-in 
technique. We used small bias voltages {U = 20 mV) and 
moderate currents (/ = 3nA) to minimize the influence 
of the tunnel process on the field-dependent magnetic 
evolution within the PdFe bilayei^^. 

PdFe/Ir(lll) bilayers have a uniaxial anisotropy of 
“ea sy-axi s” type and exhibit chiral modulations of Neel- 
typ^29l30]_ It was also established by SP-STM observa¬ 
tions that a cycloid (Fig. (b)) is the ground state of 
PdFe/Ir(lll) fllm^^. The material parameters of model 
([^ for PdFe/Ir(lll) at T = 4.2 K determined ii]SSIlII 
yield the following values for the characteristic param¬ 
eters ( p!^ : Ld = 6.44 = 23.85ao nm i.e. po = 24 
(ao = 0.27 nm is the lattice constant ), poHjj = 3.46 T, 
Ko = 1.9 X 10^ J/m^ Kd = poM‘^ /2 = 0.76 x 10^ J/m^ 
(“shape anisotropy”). The sufflciently strong values of 
“easy-axis” anisotropy {k = KjK q = 1.315) ensures the 
stability of chiral modulatio ns in PdFe/Ir(111) Aims with 
respect to stray-field effect^^^^ and make them conve - 
nient objects for investigations of chiral skyrmion d^^ l ^^ i 

The calculated magnetic phase diagram of easy-axis 
chiral ferromagnets includes the existence areas of one- 
dime nsion al modulations and skyrmion lattices (Fig. 

Inset:pIIIl. 

These chiral modulations and transitions 
between them have been directly observed by Lorentz 
transmission electron microscopy (LT EM) in fre e stand¬ 
ing nanolayers of cubic helimagnet j^^ * ^^ * ^^ * ^^ * and in 
PdFe/Ir(lll) films by SP-STIVP. 

Isolated skyrmions and their internal structure 


have been investigated by SP-STM in PdFe/Ir(lll) 
bilayerFollowing Ref.l^ we reconstruct the mag¬ 
netization profile 0{r) for one of the isolated skyrmions 
in the Aim at the applied field poH = 1.11 T (Fig. [^. 
These experimental results are in a close agreement with 
the solution of Eq. ( [I7| ) for k = 1.315 and h = 0.321 (or 
PqH = 1.11 T. In free-standing Alms of magnetically soft 
cubic helimagnets, chiral skyrmions readily con dense into 
hexagonal lattices below hs{k) (Eig. inset 
At low temperatures, however, an enhanced coercitiv- 
ity of PdEe/Ir (111) bilayers prevents the formation of 
skyrmion lattices below hs{k) (see the results of Ref.l^ 
for T = 4.2 K). This oAers a unique opportunity to inves¬ 
tigate isolated skyrmions in a broad range of the applied 
Aelds. 


Eigurej^shows selected frames from the whole SP-STM 
data set where the evolution from isolated skyrmions 
at high Aelds to spin spirals at zero Aelds can be ob¬ 
served. The two-lobe appearance of skyrmions is du e to a 
predominantly in-plane magnetization of the Cr tipP^^^. 
The strip-out of skyrmions starts in Eig. (c) where 
a skyrmion, labeled (1), has jumped to a diAerent po¬ 
sition and a skyrmion (2) has developed an elongated 
shape. In Eig. (f) more skyrmions have adopted elon¬ 
gated shapes, a process that seems to be influenced and 
assisted by defects, see skyrmion (3), and the repulsive 
interactions with other skyrmions and chiral modul ations 
along the sample edges (so called surface twistsj^^^. 
The strip-out process can be quantifled more accurately 
in an area with only one strongly pinning defect, see de¬ 
tailed view in Eig. 10 In Eig. (c) the skyrmion shape 
starts to deviate from rotational symmetry at poH 1.10 
T and becomes more and more elongated through Eig. 
(d) and (e). Other skyrmions retain axial symmetry even 
at much lower Aelds. The calculated value of the strip-out 
Aeld for k = 1.315 equals poHei = 0.65 T. The images in 
Eigs. [^(see also video materials in Ref.^^ show that the 
elliptical instability Aeld has diAerent values for diAerent 
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FIG. 10. (color online). Isolated skyrmion core gradually in¬ 
creases with decreasing field above the critical strip-out value 
Hei (a,b) and stretches into a spiral domain for H < Hei (c-f). 


skyrmions and strongly depends on skyrmion-skyrmion 
interactions, interactions with sample edges and defects. 
Similar effects are characteristic for strip-out instabilities 
of isolated bubble domains (see e.g.S^. 

To determine the skyrmion collapse field with reason¬ 
able statistical accuracy, we have monitored the repeated 
creation and annihilation of a skyrmion at higher tunnel 
bias and current as a function of applied field. With the 
tip positioned above the pinning defect, we monitored the 
telegraph noise in the spin-resolved dl/dU signal and ex¬ 
tracted the average lifetime of the skyrmion as a function 
of applied field (see insets in Fig. pT| . The skyrmion life¬ 
time decreases roughly exponentiall}!^ down to a value 
of 5 ms at 4.5 T. Minimization of functional (27) with 
Po = 24 and k = 1.315 yields the collapse field 4.4 T 
(cf. with collapse field of 7.1 Tesla calculated for zero 
anisotropy 


IV. CONCLUSIONS 

Detailed SP-STM investigations of magnetic states in 
PdFe/Ir (111) thin films and a comprehensive theoretical 
analysis within the standard model (§ enable to describe 
the basic magnetic properties of isolated chiral skyrmions 
and describe their evolution in a broad range of applied 
magnetic fields. 

The equilibrium states of isolated axisymmetric 
skyrmions are described by differential equation 0 
common for different groups of chiral magnetJ^. More¬ 
over, similar equations describe axisymmetric solitonic 
states in other co ndens ed matter systems with broken in¬ 
version symmetr}ff^J^^H^. This implies a universal char¬ 
acter of chiral skyrmion properties and allows to con¬ 
sider the investigations in PdFe/Ir films as representative 
of the entire phenomenon. These investigations include 
general features of the chiral skyrmion evolution in the 
applied magnetic fields terminated at lower fields by the 


FIG. 11. Skyrmion lifetime as a function of external magnetic 
field. SP-STM tip is held stationary above a skyrmion posi¬ 
tion while tunneling with U = +600 mV, / = 100 nA. This 
induces a continuous, stochastic switching of the magnetic 
state under the tip between skyrmion state (1) and ferromag¬ 
netic state (0) (see resulting telegraph noise in insets). Data 
points show the average lifetime of the skyrmion state as de¬ 
rived from the telegraph noise signal. At 4.5 T, the lifetime 
approaches the time resolution limit of the STM while still 
detectable. Therefore the skyrmion state is still metastable 
and the collapse field must be higher than 4.5 T. (T = 4.2 K) 


formation of skyrmion condensates or by elliptic insta¬ 
bilities of individual skyrmions and the collapse of the 
skyrmion core at high fields. 

In this paper we have investigated magnetic proper¬ 
ties of solitary skyrmions only and neglected their inter¬ 
actions with other skyrmions, with chiral modulations 
arising at the sample edge^^, and different types of de¬ 
fects. We also have considered skyrmions magnetically 
homogeneous along their axis. This assumption is justi¬ 
fied by the structure of PdFe/Ir (111) bilayers that ex¬ 
clude magnetic modulations along the film thicknesJ^. 
In thicker cubic helimagnets and uniaxial ferromagnets 
with Dn and Cn symmetry, however, such modulations 
are physically admi ssible and influence magnetic proper¬ 
ties of these system^^^^. 
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Appendix 


2. Solutions for one-dimensional modulations 


1. Skyrmion structure in different classes of 
uniaxial helimagnets 

Functional forms oiwD energy contributions are deter¬ 
mined by crystallograp hic sy mmetry of a noncentrosym- 
metric magnetic crystaP^^^. 

± ± , (A.l) 

Cnv--[w^^^], D2d--[w^]^\ Dn.[w^2\^xv\^ 

(A. 2 ) 

For e{p,ip), 'tp{p,(p) 

= cos{'ip — ‘p) 0 p — sin 6* cos 0 sin 
— sin ('0 — ip)0^ H— sin 0 cos 0 cos('0 — (f)'ipcp (A.3) 


For one-dimensional modulations propagating along 
the axis the energy functional (© can be written as 

n; = A( 6 >|+sin^ 02p‘^)^WD—Kcos^O — /j.oMHcosO. (A.9) 


In the DM energy contribution wjj (A.9) Lifshitz in¬ 
variants of type induce modulations propagat¬ 

ing in the xy plane (e.g. helicoids and cycloids in Fig. 

(z) 

2 (a), (b)), and invariants Cxy favour modulations along 
the 2 ;-axis {cones). 

Helicoids and cycloids. To be specific, we consider 
in-plane modulations propagating along the x— axis. 
Depending on the magnetic crystal symmetry, different 
types of modulations are stabilized by the wd energy 
functionaP. Particularly, in cubic helimagnets and uni¬ 
axial ferromagnets of Dn crystallographic classes, M ro¬ 
tates as a Bloch-type domain wall {helicoids)^ and in uni¬ 
axial ferromagnets with Cnv symmetry the magnetization 
rotates along the propagation direction like a Neel-type 
domain wall {cycloids) (Fig. 1 (a), (b)): 


^ = cos ('0 + ^) 0 p — sin 0 cos 0 sin ('0 -h 
— sin('0 -h ^)0^p -sin 0 cos 0 cos{'ip -h (A.4) 


m = Cy sin^(x) -h cos^(x) (helicoids), 

m = Cy sin^(x) -h cos^(x) (cycloids). (A. 10 ) 

The Euler equation for the functional 


= sin('0 — (f)0p — sin ^ cos ^ cos (V^ — (p)ipp 

+ - cos('0 — (p)0u^ -sin ^ cos ^ sin('0 — ^)'ip^ (A.5) 

P P 

W 2 ^ = sin('0 -h p)0p — sin 0 cos 0 cos{^ -h ^)'0p 
+ - cos('0 + p)0y^ — - sin6>cos6>cos('0 + p)^^^ (A.6) 

The solutions '0(0) are determined by crystal classes 
of the systenP. 

Cnv : 0 = D2d : 0 = -(^ + 7r/2, Dn ' = p + 7r/2, 

5'4 : 0 =-(^-h0i, Cn : 0 = (/:? +01. (A.7) 


WY,{0) = AOl - DOx - poMH cos 0 - K cos^ 0 (A. 11 ) 


yields magnetization profiles 0 {x) for helicoids and 
cycloidP. Analytical solutions for 0{x) describe heli¬ 
cal modulations distorted by the applied field and uni¬ 
axial anisotrop 3 P. These helicoids (cycloids) gradually 
unwin d into a set of isolated domain walls at critical line 

Cones. In cubic helimagnets and uniaxial ferromagnets 

with Cn and Dn symmetries, the DM energy functional 

(z) r 

Wd includes Lifshitz invariants Cxy (A. 2 , A. 8 ) favour¬ 
ing chiral modulations {cone phases) along the z— axis. 
Minimization of the energy functional 


wi^{0) = AOl — DOz — PoMHcosO — K cos^ 0 (A. 12 ) 


For ferromagnets belonging to S 4 and Cn classes en¬ 
ergy functionals wd include two terms: wd = Diw[~^ -f 

D 2 W 2 ~^ for S 4 and wd = Diw[~^^ + D 2 W 2 ^^ for Cn- An¬ 
gles 01 = arctan(D 2 /Di) and the effective values of the 
DM constant are D = D‘1 + 

For noncentrosymmetric cubic ferromagnets belonging 
to T and O crystallographic classes the energy functional 
Wd has the following fornP^ 

Wd = + Ipy} + = m ■ rotm, (A.8) 

and stabilizes solutions with 0 = 7 r /2 + (/:?. 

The skyrmion energy densities for all these structures 
can be reduced to a common functional fornP. 


yields the solutions for singl e-harmonic modulations de¬ 
scribing the cone phastPESl 


cosP = 


H 



0 = z/Ld- (A.13) 


These equations include the characteristic parameters of 
a uniaxial chiral ferromagnet (16). 


3. Characteristic lengths and critical parameters 

In uniaxial noncentrosymmetric ferromagnets, chiral 
modulations arise as a result of a competition between 
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FIG. 12. Projections of m onto the basal plane in the core of chiral skyrmions for noncentrosymmetric uniaxial ferromagnetfl 


the DM interactions favouring a rotation of the magne¬ 
tization, the exchange coupling and the “potential” en¬ 
ergy f{0) = —fioMH cos 0 — K cos^ 0 tending to sup¬ 
press such modulations. The balance between the chi¬ 
ral energy wniO^Op) and the potential energy contribu¬ 
tions f{0) determine the equilibrium spin configurations 
in chiral magnets. At zero field and for zero anisotropy 
single harmonic modulat ions 0 = qqx (go = D/{2A)) 
minimize the functional ( |A.9| ). In the opposite limit of 
strong anisotropy [K > 7r^D^/(16A)) these modulations 
transform into a set of isolated 180° domain walls sepa¬ 
rating the homogeneous states with 0 i = 0 , O 2 = tt. 

The width of an isolated Bloch domain wall its 
energy, 7 ^, and anisotropy field Ha are as follows: 

= (AA4) 

These are the fundamental parameters describing mag¬ 
netic states in a common (centrosymmetric) uniaxial 
ferromagnefP^. To demonstrate a competing character 
of the magnetic interactions in chiral uniaxial ferromag- 
nets, we consider an isolated domain wall at zero field 
that separates the homogeneous states with 0 i = 0 and 
$2=71. The equilibrium states of thi s dom ain wall are 
derived by minimization of functional ( |A.9| ) for H = CP^. 
The standard calculation of the wall energ:\P ^ 7 ^ = 
/o^ ['^h(^) — '^h(O)] dx yields the following resulfP^ 


7 w = 4yA^-7r|D| .(A.15) 


The first (positive) term in (A.15) is the wall energy of 
a uniaxial ferromagnefP^ arising as a common effect of 
the uniaxial anisotropy pinning the magnetization vector 
along the easy-axis and the exchange stiffness suppress¬ 
ing deviations of M from these directions. The negative 


energy contribution in 7 ^ (A.15) is due to the DM inter¬ 


actions favouring modulations of the magnetization with 
a specific rotation sense. The strength of this “wind¬ 
ing force” is characterized by XjLjj'. the larger the DM 
coupling, the smaller the period of the modulations. For 
L]j < ttLb the wall energy becomes negative manifesting 
the instability of the homogeneous states with respect to 
chiral modulations. 

The dimensionless parameter x introduced 


7v 


= 4 \/AK (1 — x), 


TT 

K = — 


\D\ 


'kLi 


4 x/AA Lb 


, (A.16) 


provides the criterion for the existence of chiral modu¬ 
lated states. For x > 1, the DM interactions overcome 
a pinning of the magnetization along easy-axis direction 
and stabilize modulated states. For 0 < x < 1 chiral 
modulated phases are totally suppressed, and chiral pat¬ 
terns exist as metastable localized states in a form of 
isolated skyrmions and domain walls (kinks). 

Parameter k ( |A.16[ ) is similar to the Ginzburg-Landau 
parameter kql = \/i in the theory of Abrikosov vortices 
(the mixed state) in superconductor^^. The parameter 
>cgl is a ratio of two characteristic lengths, the eoherenee 
length ^ and the penetration depth A. Abrikosov vortices 
exist in super condu ctors with hql > l/\/2 {type-II su- 
pereonduetors'^^^. Physical analogies between super¬ 
conductor’s mixed s tates and chiral magnetic modula¬ 
tions are discussed irP^. The characteristic lengths Lb 
and Lb provide different ways to introduce reduced vari¬ 
ables into model (12). 


The reduced energy functional based on the charac¬ 
teristic length Lb with control parameters h = HjHa 
and K (A.16) a nd th e reduced magnetic field (where Ha 
define in Eq. ( |A.14 ) is convenient for investigations of 
isolated skyrmions in helima gnets with a strong uniaxial 
anisotropy (0 < x < 
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